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We include effects of nonlocal quark condensates into QCD sum rules (QSR) via the Kallen- 
Lehmann representation for a dressed fermion propagator, in which a negative spectral density 
function manifests their nonperturbative nature. Applying our formalism to the pion form factor 
as an example, QSR results are in good agreement with data for momentum transfer squared up to 
fa 10 GeV^. It is observed that the nonlocal quark condensate contribution descends like 1/Q^, 
different from the exponential decrease in obtained in the literature, and contrary to the linear 
rise in the local-condensate approximation. 

o 

^ ■ PACS numbers: 11.55.Hx, 12.38.Aw, 12.38.Bx, 13.40Gp 

^ ■ In QCD sum rules (QSR) nonperturbative contributions are taken into account via vacuum expectation values of 
nonlocal operators, such as {q{0)q{z)) and (G(O)G'(z)) where g is a quark field and G is the gluon field strength. 
In the standard approach vacuum effects are assumed to be sufficiently soft to allow the Taylor expansion of, for 
instance, the quark condensate (g(O)q(z)), at z = by means of local composite operators, 

^ A local condensate (qq), i.e., the first term of the above expansion, prohibits momentum flow. A loop diagram then 
turns into a tree diagram as shown in Fig.[T] when inserting the local quark condensate into the lower (nonperturbative) 
line. The external momentum q flows only through the upper (perturbative) line, and one has the loop integral 
CN ■ approximated by the product of the propagator 1 /q^ and the condensate (qq) [2] . With this localization assumption, 
^ ' simple hadronic properties including masses, decay constants, moments of hadronic wave functions, and form factors 
L have been calculated in QSR. 
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FIG. 1: Loop diagram with the insertion of the local quark condensate, where k denotes the loop momentum. 

It has been known that nonperturbative contributions from local quark condensates grow with the momentum 
transfer squared in form factor calculations, whereas perturbative contributions decrease [1, 1^. This is the 
reason why the standard QSR approach encounters difhculty, when applied to form factors in the region with high 
> 3 GeV^ ^5|]. It has been observed that the dependence of nonperturbative contributions is moderated by 
employing the nonlocal quark condensate {q{Q)q{z)) Moreover, using local quark condensates in QSR analysis of 
more complicated processes, such as Compton scattering which involves four-point correlation functions infrared 
divergences appear. Consider the box diagram in Fig. [21 where a light hadron is scattered by an on-shell photon 
of momentum qi. The external momentum qi flows through the upper horizontal quark line, when the local quark 
condensate is inserted into the left vertical quark line. The upper line then gives a divergent propagator proportional 
to l/ql — )> cxo, and the evaluation of the Wilson coefficient associated with the quark condensate makes no sense. A 
resolution of the above difficulties is to relax the localization assumption. Including the nonlocal condensates, a finite 
loop momentum k is allowed to fiow through the box diagram, and the above infrared divergence is smeared into 



* E-mail: hsiehrcCSphys. sinica. edu.tw 
t E-mail: hnliCSphys. sinica. edu.tw 



l/{qi + fc)^. This is our motivation to investigate effects of the nonlocal quark condensates in QSR. In this letter we 
shall set up the framework by studying simpler processes like the pion form factor, and compare the results with the 
local condensates and with the nonlocal condensates. 
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FIG. 2: Box diagrams for Compton scattering. 

Nonlocal condensate models Q have been applied to QSR for the pion wave function Isj-fllj. whose outcome was 
then treated as an input of the perturbative QCD factorization formula for the pion form factor Recently, 
Bakulev, Pimikov and Stefanis calculated the space-like pion form factor based on QSR with nonlocal condensates 
0. They parameterized the nonlocal quark condensate as (g(O)q(z)) = {qq) exp{—\z^\Xg/8) where is related 
to the average virtuality of the condensed quarks. Our formalism is different, which starts from the Kallen-Lehmann 
(KL) representation for a dressed propagator of the quark q (l^ . 
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where represents the exact QCD vacuum, T denotes the time ordering, the spectral density functions 2(/^^) 
describe the glutinous medium effect, and /i is the effective mass. The KL representation can be deemed as a 
superposition of free quark propagators for all mass eigenstates with the weights pf 2(p'^)- 
Equation ([2]) is recast into 



{n\T{q{z)q{Om) 
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ds exp [ — s 
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(3) 



We decompose the above matrix element into the perturbative and nonperturbative pieces 



{n\T{q{z)q{0)m) = iZS{z,mg) + {n\ : q{z)q{0) : |r!), 



(4) 



respectively, with Z being a renormalization constant, S{z, ruq) being the quark propagator in perturbation theory, 
niq being the quark mass. The nonperturbative piece collects the contribution from large /i^. 
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The lower bound for the integration variable p^ is usually set to the multi-particle threshold in the KL represen- 
tation. In the above expression we have modified it into pj. — max(cs, m^), where the free parameter c of order unity 
will be fixed later. This choice respects the multi-particle threshold, and at the same time guarantees a finite integral 
in Eq. Note that the integration over /x^ in Eq. ^ develops a divergence as the variable s approaches infinity 
without the above modification. A negative spectral density function implies confinement 15|, and we indeed have 
the property p\{p'^) < as shown in our formalism. 
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We define the distribution functions 
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and assume tlie spectral density functions to be 



dfi^ exp 
d/i^ exp 



s 
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with tlic free parameters a, Ni and A^2- It is then easy to see that 



fs{s) cx 



■exp(-^^/s - afil), 



(6) 
(7) 

(8) 
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exhibits the limiting behaviors exp(— rn^/s) at small s and exp(— acs) at large s, consistent with exp(— to^/s) and the 
exponential ansatz exp(— cr^s) proposed in the literature fill [l6j. respectively. Hence, we set the threshold mass to 
a constituent quark mass ~ 0.36 GeV [l3|. Comparing the Taylor expansion of the nonlocal quark condensates 
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mhM^)) = -Tr[7^(r!| : q{z)q{0) : \n)] 



(10) 



with Eq. ([21), we obtain the constraints 

poo 

/ fy{s)ds = mq, 
Jo 

/ fs{s)ds=l, 
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sfs{s)ds = - {Xl - ml), 



(11) 



which determine a, Ni and N2 in Eq. (|S]), given values of and niq. 

The dressed propagator includes both the perturbative and nonperturbative contributions, 
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with the definitions 
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The second and third terms on the right hand side of Eq. (112]) arise from the background gluon field [19|, |2C| , and the 
forth term comes from the nonlocal quark condensates with the integrations over p? and s being exchanged in Eq. ([SJ . 
As stated before, local quark condensates lead to contributions linear in Q^, which are more serious than the constant 
contributions from local gluon condensates at large [1, 01 • Gluon condensate contributions to the pion form factor 
are actually negligible, and quark-gluon-antiquark condensates vanish in the localization limit. Besides, the potential 
infrared divergences in four-point correlation functions are generated from local quark condensates. Therefore, only 
the nonlocal quark condensates are taken into account here. 

Inserting Eq. (1121) into the triangle diagrams for the three-point correlation function, we derive the perturbative 
and nonperturbative contributions to the pion form factor i<V((5^), 



/2F^(Q2)exp (^-^) = ^ 1^'" dsids2P^'=''*(si,S2,g')exp 



■Sl + 52 
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FIG. 3: (a) Perturbative contribution, (b) two-quark condensate contribution, (c) four-quark condensate contribution, and (d) 
gluon condensate contribution to tiie pion form factor. 



In the above expression /^r is the pion decay constant, is the pion mass, M is the Borel mass, and sq is the 
duality interval The calculation of the spectral function fp'^^^ associated with the perturbative contribution, and of 
the quark (gluon) condensate contribution A''""'"'^ (A^'""") involves four types of diagrams displayed in Fig. [3l The 
perturbative spectral function and the gluon condensate contribution are given by 0,13] 

= ^^A^-j^Q^{s^[Q^+S,f + S2{Q^+S2f-S^S2[2Q^ + Q\s^+S2)-2{s\+sl)+&S,S2]]^ (15) 



respectively, with Nc being the number of colors, Cu (e^) being the charge of the u (d) quark, and the variable 

A = (si+S2 + Q')'-4siS2. (17) 
We compute the quark condensate contribution, obtaining 

^,uark ^ ^-^^ ^^^^^^ J„ _ ^Jf)pa^^ + ^^Jf _ e./^p^^J CXp (-^^) + (g?)' A^^^' ■ (18) 

The two-quark condensate spectral functions P2gt the four-quark condensate function A|^"'^'^ are written as 

a ^ 

2qc yjl2 



-3Q^(Q6 + 2QV' - + si) - Q^iQ"^ + 2QV' + 40'/ + 3/i*^)(si + S2) 

-2(Q4 + 2gV' - 2/)(s?S2 + slsi) + g2(Q4 ^ 2Q2^2 _ 5^4)(^2^^ ^ ^2^^) 

+ 6(Q'* -f OV' - /)s?S2 - 20^(0^ + 50V' - 20V* - 6m')siS2] , 



+3(g4 - 2gV' + 2/)(s? + s^) + Q^iQ" + GQ^fJ-^ - 18//) (si + S2) 
+2{sls2 + slsi) - (Q2 _ 6fi^){sls2 + S2S1) - 6s?S2 - 2(0^ - 12gV + 6/)siS2] , (19) 

Air' = - ed/2') Jm^ I / , '^'^ J + ds, ds2 j (e-^^/*^' - l)gie-^^/^' 

+ [eJi - ej^) Jm^ ^ f ds^ ^ ds2 + £ ds^ j'^ ds\ (g-^^/*^' - 1)526^^^/^^' , (20) 



pI,c = -^[2gV(3O^-12gV + 10A*') + (4 + 4) + 3(Q'-2/.2)(.? + si) 
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with the functions 



91 = ^^^{Q'[-6m^^^^^Q^+^^}+n^^{4Q^ + 6^^^)] + {Q^ + ^^^)sl + {Q'-^,')sl 

52 = T;^?^ {Q'[Q' + - 6/i2(Q2 _ /,2) ^ 2m2(Q2 _ 3^2)] ^ (^4 _ 4^2g2 ^ q4^^2 ^ (q2 ^ ^2)2 

+2Q2s2(Q2 + m2)(Q2 + rn^ - 3^2^ _ 2Q2g^[2m4 _ (g2 _ ^2^)(g2 _ 3^2-)^ 
-2siS2(m*-m2Q2_2Q-*)}, 



(21) 



and the variables 



a = {rri^Q^ + /i^S2) 

p = (M2g2 + „j2^^^ 



1 



1 



1 1 



+ si — fi'^ 



(22) 



Note that the singularity from si (s2 0) in the function gi (172) is removed by the factor (e — 1) 

[{e~''^^^^ — 1)] in Eq. ((20)) . It is observed that the contributions from the nonlocal quark condensates must be 



power-like in in the asymptotic limit, no matter how to parameterize p\ 2(/^^)- The dominant contribution 



^4qc 



descends like l/Q^ as — > cxd, which is different from the exponential decrease in obtained in [5J, and contrary 
to the linear rise in the local condensate approximation [3!, [1| . 

The local condensates appearing in Eqs. (|16p and ([T5)) are taken to be [2l| 



TT 



0.005 ±0.004 GeV^ 

-(1.65 ±0.15) X 10^2 GeV^ 
(1.5 ±0.2) X 10"* GeV^ 



(23) 



The duality interval So(Q^) at a given is determined by the requirement that the form factor is least sensitive to 
the Borel mass M. The average virtuality Xq and lower bound c, being not known with certainty, are fixed by fits to 
the data of the pion form factor F^{Q^) = 0.179 ± 0.021 at = 1.99 GeV^ [H-IH- In figure IDJa) we display the 
allowed values of c and Xq as a curve in the c-Xq plane. The range of Xq is consistent with Xq = 0.63 GeV from QSR 
[26] and Xq = 0.83 GeV from the instanton analysis [27]. Below we adopt Xq = 0.745 GeV and c = 0.25 to produce 
the central values of our predictions for the pion form factor. Choosing the light quark masses m„ — 4.2 MeV and 
nid — 7.5 MeV, we solve for the free parameters a, Ni and N2 from the constraints in Eq. ([TT|). whose results are 
listed in Table HI The product ac w 7 GeV~^ is in agreement with the value of cTq « 10 GeV~^ postulated in [2^. 
The opposit e sig ns of A^i and iV2 imply the violation of positivity, which can be interpreted as a manifestation of 
confinement 1151. 





A, (GeV) 


rUg (MeV) 


a (GeV^^) 


N,/{qq) (GeV-*) 


N2/{qq) (GeV-*) 


u quark 


0.745 


4.2 


28 


40.32 


-18790.09 


d quark 


0.745 


7.5 


28 


72.00 


-18790.09 



TABLE I: Parameters associated with the quarks u and d in our formalism. 



Figure lU^b) indicates the best choice of sq — 0.715 GeV^ with Xq = 0.745 GeV and c = 0.25, at which the pion 
form factor Ft^{Q'^ = 1.99 GeV^) becomes independent of M for M > 1.5 GeV. In the calculation below, we simply 
set the Borel mass to M — 1.5 GeV. In Fig. |3^c) we present the dependence of the best choice So(Q^) with the 
same inputs, whose curve is close to a straight line: 

so(q2) ^0.6 + 0.06Q2_o.002Q^ (24) 

for > 1 GeV^. It is seen that Sq drops rapidly in the region of low < 1 GeV^, where QSR are supposed to 
be inapphcable. sq in Fig. I^Jc), increasing from 0.65 GeV^ to 1.05 GeV^ for 1 GeV^ < < 10 GeV^, shows a bit 
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stronger dependence compared to that in [5|. Nevertheless, its range obeys the postulation Q that it should not 
be lower than the middle point 0.6 GeV^ of the interval between the meson masses = and m^^ = 1.6 GeV^. 
Besides, we have confirmed that the pion decay constant squared takes the value ~ 0.0171 GeV^ for sq ~ 0.7 GeV^ 
in our formalism with the nonlocal quark condensates. 




M(GeV) Q=[GeV=] 



FIG. 4: (a) Curve for the allowed values of c and \q that produce the data of Ft,{Q^ = 1.99 GeV'^). (b) M dependence of 
F^{Q'^ = 1.99 GeV^) for different so with A, = 0.745 GeV and c = 0.25. (c) dependence of so for the pion form factor. The 
function of the fitting curve is presented in Eq. (|24|) . 

Our results of the pion form factor F-^{Q'^) are displayed in Fig. EJa) for three values of \q = 0.765, 0.745, and 
0.725 GeV, corresponding to the curves from top to bottom, respectively. Their difference indicates the theoretical 
uncertainty of our analysis. It is obvious that all three curves are well consistent with the experimental data for 
> I GeV^, the region where QSR are applicable. We investigate the perturbative and condensate contributions 
to the pion form factor F^iQ"^ = 1.99 GeV^) at different Borel mass M, as shown in Fig. [5jb). It is observed that 
the former increases with M, and the latter decreases with M for M > 1 GeV. The gluon condensate contribution 
becomes negligible for M > 1 GeV, justifying the sole modification from the nonlocal quark condensates. Although 
the magnitudes of different contributions vary with M, their sum is almost constant for M > 1.5 GeV. The quark 
condensates contribute 23% of the pion form factor F^{Q'^ = 1.99 GeV^) at M = 1.5 GeV, which is slightly higher 
than the percentage 17% in the localization approximation Q. 

In this letter we have included the nonlocal quark condensates into QSR via the KL parametrization for a dressed 
fermion propagator, which is decomposed into the perturbative and nonperturbative pieces. The negative spectral 
density function implies that the contribution from higher effective quark masses is nonperturbative. The parametriza- 
tion of the spectral density functions leads to the known exponential ansatz for the nonlocal condensate model in our 
formalism. We have analyzed the pion form factor as an example, and the results are in good agreement with the 
data for between 1-10 GeV^. The fitted virtuality = 0.745 GeV and the derived duality interval so(Q^) are 
also consistent with those reported previously. The nonlocal quark condensates remedy the improper dependence of 
the nonperturbative contribution in the localization approximation at large Q^: the quark condensate effects decrease 
like l/Q^, which is different from the exponential decrease obtained in the literature. Viewing the success of this 
approach to the pion form factor, we shall extend it to more complicated processes, including Compton scattering 
and two-photon hadron production [29j . 
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in July, 2008. This work was supported by the National Center for Theoretical Sciences and National Science Council 
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FIG. 5: (a) dependence of for, from top to bottom, Aq = 0.765, 0.745, and 0.725 GeV. The data points are referred to 
[l^-dH]. (b) M dependence of the perturbative and condensate contributions to Ftt{Q^ = 1.99 GeV^). 
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